Simulation of Current Experiments on the Tunneling Effect 
of Narrow e-Near-Zero Channels 
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In this paper, we discussed the current experiments on the tunneling effect of electromagnetic 
energy through narrow channels of e-near-zero(ENZ) medium in the microwave range. Using the 
finite element method, we carried out a full wave simulation of the two kinds of experimental 
configurations at present. It was shown that the ability of the electromagnetic waves penetrating 
into the ENZ medium is very necessary. The present experimental setups only using metamaterial- 
filled narrow channel in the configuration of parallel plate waveguides are unlikely to realize effective 
tunneling. Contrarily, the artificial plasma medium emulated using hollow metallic waveguide can 
achieve nearly perfect tunneling in a narrow channel without any transition section around its cutoff 
frequency, which exceed the scope of original ENZ tunneling theory and can be described by a simple 
equivalent circuit model. 

PACS numbers: 78.20.Ci, 41.20.Jb, 52.40.Fd,78.66.Sq, 42.82.Et 



I. INTRODUCTION 

Since the rediscovery of the theory first proposed by 
Veselago in 1968 0, metamaterials have been an attract- 
ing subject of growing worldwide interest. The metama- 
terials can possess many exotic electromagnetic proper- 
ties which are absent in natural materials. So they are of 
particular importance in the engineering of electromag- 
netic wave fields, which is leading to a variety of new 
microwave or optical devices. 

Most current research is still focused on the double neg- 
ative (DNG) metamaterials in which the permittivity and 
permeability can both less than zero simultaneous, lead- 
ing to the negative index of refraction. Recently, there 
is a surge of investigation on e-near-zero(ENZ) medium 
whose effective permittivity e is near zero and effective 
permeability fi remains natural. In 2006, Silveirinha and 
Engheta0 have theoretically predicted that the electro- 
magnetic waves can be squeezed and tunnel through very 
narrow subwavelength channels filled with ENZ medium, 
independent of the specific geometry of the channel. The 
incoming wave front can be replicated at the output in- 
terface. The reflection at bends or junctions of different 
devices can decrease greatly and even disappear using 
this method. So, the ENZ medium will play an important 
role in the applications of transport and interconnects of 
microwave and optical wave energy. 

Till now, several experiments have been reported to 
testify Silveirinha and Engheta's ENZ theory and the 
seemingly 'counterintuitive' tunneling effects are demon- 
strated at the certain microwave frequencies [3, H, H, 0. 
The tunneling geometries in these experiments are ba- 
sically 'U-shaped'. At present, there are two ways to 
experimentally realize the ideal ENZ tunneling channels 
proposed by Silveirinha and Enghcta. 

In the work by Liu and Cheng et al0, the rela- 
tively complex artificial resonant inclusions implanted in 
the parallel plate waveguides or microstrip circuits are 
adopted to achieve the effective ENZ medium near the 



resonance frequency. The two-dimensional(2D) scatter- 
ing of the electromagnetic wave in the original theory 
are simulated by the incident polarized nearly transverse 
electromagnetic(TEM) waves. 

In the later experiments, the other way was adopted 
by Edwards et al [H, 0. The artificial plasma medium 
emulated by the TEio guided mode propagating in hollow 
rectangular metallic waveguides, mimics the equivalent 
response of ENZ materials near the cut-off frequency of 
guide wave modes. 

In all above experiments the transition section in which 
the incident mode "penetrates" into the narrow channel 
is less considered. Intuitively, the sharp discontinuities 
will cause strong reflection of electromagnetic wave and 
excite high order evanescent modes, which not consid- 
ered in the discussion in current work. To investigate this 
problem, we performed the full wave simulations of these 
experiments respectively in the following sections. From 
these results, we emphasis the importance role of the 
transition section in the tunneling effect of the 2D ENZ 
channel. The experiments results must be cautiously 
examined to exclude other possibilities before being af- 
firmed to the ENZ tunneling phenomena predicted by the 
theory. While for the case of hollow metallic waveguides, 
The effective tunneling can be realized and the transi- 
tion section is unnecessary. The corresponding pictures 
of underlying physics are discussed. 



II. ENZ CHANNEL BETWEEN PARALLEL 
PLATE WAVEGUIDES 

Over the past years, the development of metamaterials 
has greatly expanded the realm of available permittivity 
e and permeability \x [|[ . The artificially microstructured 
composites, may have become the most straight way to 
realize the effective medium with certain electromagnetic 
roperties, including ENZ. In the work by Liu and Cheng 
(, the planar structure of complementary split ring 
resonator (CSRR) are adopted, which have an electric 
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FIG. 1: (a) A U-shaped geometry of the microwave tunneling 
problem in a 2D parallel plate waveguide theoretically dis- 
cussed by Silveirinha and Engheta in Ref[3[. (b) Setup of our 
simulation, PMC boundaries conditions are on the x=0 and 
w planes. Other outer faces are assigned as PEC boundaries. 
a = w = 11mm, and a c i=lmm. (c) Simulation results of (b): 
the transmission coefficients |jSai| of ENZ channels with the 
same geometry (Lo = 15mm) and different permeabilities 
(jj./ fj-o = 1) 2). The one for a air-filled channel is also plotted. 
As a comparison, l&il for the ENZ(/i//io = 1,2) tunneling 
model of the same size with side ENZ transition section 
[L\ = Z/2 = a c h) is plotted as dashed lines. 



resonance in the certain range of frequency [9]. The res- 
onance permittivity values of effective medium can be 
achieved. 

A general U-type geometry for the 2D tunneling prob- 
lem of a ENZ channcl(shown in Figure 1(a)) has been 
discussed in detail as a prototype in detail in Silveirinha 
and Engheta's theory Q. Although contributing to the 
total cross section (A p — (L±ai + L^a^) + La c h) and in- 
conveniencing the approximate condition (ko{j, r A p — ► 
in Ref0)to realize effective tunneling, the lateral ENZ 
sides is retained as transition section, allowing the elec- 
tromagnetic wave to penetrate into the narrow channel. 

But in current experimental setup [3, Q, only chan- 
nel filled with metamaterials are considered, without any 
transition section of effective ENZ medium. The lateral 



metallic step sides are between parallel plate waveguides. 
Although in the paper by Liu et al, they are considered 
as a shunt admittance in a simple equivalent model [3]. 
But in their approximation, it was oversimplified and 
can only be valid in a very limited frequency range. 
The inhomogeneities at the metallic step in waveguides 
are much more complicated. Furthermore, the incident 
TEM-mode waves in parallel plate waveguide can excite 
more high order evanesce modes, at each sharp metallic 
step, leading to the interaction between these disconti- 
nuities, which can't be reckoned in such a simple model. 
As we know, whether the tunneling effect can be real- 
ized in this kind of configuration has not been discussed 
theoretically yet. 

Now we first consider the simplified ideal model, as 
shown in Fig. 1(b). Here the two planes, x = and 
w, are assigned as perfect magnetic conducting(PMC) 
boundaries. Other outer faces except two waveports are 
assigned as perfect electric conducting(PEC) boundaries. 
These constraint conditions insure the incident waves are 
of transverse electromagnetic(TEM) mode with certain 
polarization as marked in (b) . The field invariance along 
the z direction can be fulfilled, so it degenerates to a 
2D problem independent of the width, w. If the effec- 
tive medium theory holds well for the volume with the 
CSRR or other composite structures, this will be a ideal 
approximation for the experiments by Liu and Cheng et 
al0,l. 

In our full wave simulation, the fixed height of the ENZ 
channel a c / l (=0.5mm) is much less than the y spacing 
of the input and output parallel plate waveguides, a = 
11mm. Simulations of this parallel waveguides geometry 
were performed with Ansoft's High Frequency Structure 
Simulator (HFSS) based on finite element method(l8|. 
Two set of constitutive parameters (/j, r =l, 2) for lossless 
ENZ materials and are used. 

For a fixed channel with LO = 15mm, the simulated 
results of the transmission coefficients IS21I are shown in 
Fig. 1(c) (the lines with symbols), the air- filled channel 
without lateral are simulated, and corresponding IS21I 
also plotted for comparison. 

The IS21 1 curves of the two ENZ channels without tran- 
sition section are far below unity in the most range of 
microwave frequency, though they are above the level 
of air-filled channel. For comparison, we simulated the 
ENZ tunneling model shown in Fig. 1(a) with side 
transition ENZ sections (L\ = L2 = a c h)- The IS21I 
curves (dashed lines) suggest nearly perfect tunneling 
(I £21 1 ~~ ¥ 1) through the the narrow channel with the 
same size, independent of the permeabilities. 

In the Fig. 1(b), the tangential components of electric 
fields must be zero at the vertical surfaces of both side 
steps, which is strongly non-compatible with the propa- 
gating TEM mode in the two parallel plate waveguides. 
So the discontinuities obstacles the propagation of elec- 
tromagnetic waves for air or ENZ-filled channel. On the 
contrary, the normal components of electric fields must 
be zero due to finite electric flux density inside the ENZ 
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surface, while the tangential ones are continuous across 
the interface. Thereby the lateral ENZ side layer let the 
normally incident TEM waves penetrating into the ENZ 
directly with no need of changing its configuration as 
shown in Fig. 1(a). 

Another prominent character of the IS21I curves in the 
two ENZ channels without transition sections is the IS21I 
became very significant (near unity) around certain fre- 
quencies. The position of transmission peaks in ENZ 
channels depends strongly on the constitutive parame- 
ters, 12.8GHz i.e. andl0.6GHz for the ENZ medium 
with u—1 and 2 respectively. The transmission peaks 
also exist in air-filled channels at different frequencies, 
and These have studied in detail as the Fabry-Perot 
oscillation [ijj]. Does the peaks of IS21I in the ENZ chan- 
nels have the same origin? We further simulated six cases 
with the different channels lengths, Lq = 10to60mm with 
the same permeability (/i=2). And the results are shown 
in Figure 2(a). 

According the theory of Silveirinha and Engheta, the 
anomalous tunneling effect is independent of geometry if 
only the ENZ channel is narrow enough. Here,the fre- 
quency of these transmission peaks increases as we de- 
crease the length of channel section, which shows a strong 
dependence of geometry. (The resonance peaks also alter 
as we change the width of ENZ channel, a c h (results are 
shown here) . So the origin of the tunneling effect can be 
excluded. 

To survey the underlying propagation processes in 
ENZ channels at the resonance frequencies, the distribu- 
tion of the calculated electric field at phase are depicted 
in Figure 2(b) for a typical channel with LO = 15 at the 
corresponding resonance frequency (10.6GHz). The mag- 
nitude of the calculated electric field along the central 
line of the channelled line in Figure 1(b)) is also plotted 
in the Fig 2(c). The Fabry-Perot-like character of the 
strong "standing" waves is observed in the ENZ chan- 
nel with the maximum magnitude at the two open ends 
of the channel, very similar to Fabry-Perot oscillation in 
the air-filled slots or channel [l3l|. 

Commonly, the mth-order Fabry-Perot resonance oc- 
curs at the frequencies near f =mc/(2nL0)(here c is the 
light velocity in vacuum, h is the length) as the channel 
width is small enough. The transmission peaks will 'blue 
shift' as the length channel increases, which is consistent 
with our results. So these transmission peaks could be 
attributed to the Fabry-Perot-like resonance. 

There is a non-zero minimum of the magnitude of field 
in the middle of the channel, suggesting a slight differ- 
ence from rigid standing wave picture of classical Fabry- 
Perot resonance. It may seem more confusing because 
the refraction index (n = ^fefi) of a ENZ medium is near 
zero, and the wavelength of electromagnetic waves inside 
ENZ medium is 'infinite'. So no obvious phase difference 
should exist in ENZ medium. The further thorough in- 
vestigation on these anomaly resonance peak is being on 
the way. 

These geometrical resonance characteristic of the 
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FIG. 2: (a) Simulated results of the transmission coefficients 
j^il of ENZ tunneling channels with different length Lo 
at a fixed phase (b) Calculated field distribution of electric 
vector for the Lo = 15mm ENZ channel at the corresponding 
resonance frequency, 10.6GHz. (snapshot in time). The 
corresponding magnitude of the electric field along the 
central line of the channel are shown in (c) . 



transmission peak is conceptually totally different from 
the tunneling effect associated with ENZ channels which 
is independent of specific geometry [2j. So in practice 
researchers can easily differentiate the resonance trans- 
mission peaks due to either the ENZ tunneling effect or 
these Fabry-Perot-like one by altering the experimental 
geometry, such as the channel length Lo- 
in these present experiments 0, H|, the inhomo- 
geneities also exists at the metallic vertical steps and 
the observed transmission IS21I at the passband is far 
below the perfect tunneling level. The geometries of the 
tunnel channel were fixed in their experiments. So the 
observed passband may not be identified as the results 
of ENZ tunneling. The geometries of the tunnel channel 
were fixed and never changed during their experiments. 
Thereby their results may not be the strong evidence of 
ENZ tunneling as they declared. 

A full wave ab initio simulation is necessary to give 
a more clear picture. The configuration in our simula- 
tion (shown to the left of Figure 3(a).) are nearly the 
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FIG. 3: (a) Retrieved effective permittivity for a CSRR unit 
cell(inset: the geometry of a CSRR cell). The simulation con- 
figuration for the tunneling experiment in RefQ are shown 
to the left of (a). The corresponding simulation results of 
transmission coefficients IS21I as a function of frequency are 
depicted in (b) and (c). Each tunnel channel consist of m x 
n uniform CSRR cells (m cells in the propagation direction 
and n cells in the transverse direction). The IS21 for control 
sample is also plotted for comparison. 



The un-patterned copper-clad FR4 board was used in the 
'control' channel with the same geometry, which was also 
simulated and corresponding |S 2 i| plotted for compari- 
son. Apparently, the calculated IS21I curve for channel 
with CSRR patterns is similar to the control one. But 
the resonance peak shifts about 2GHz to lower frequency 
compared to the unpatterned control sample. There are 
weak fine resonance structures between 6 and 8 GHz, 
is consistent with the electric resonance found in single 
CSRR cell retrieval procedure (shown in Fig. 3(a)). 

To find whether there exists strong tunneling evidence 
in the CSRR channel, we examined CSRR channels with 
different length(L )- Three cases, 3x4, 5x4, 7x4, are 
simulated, in which 3,5 and 7 cells are along the longi- 
tudinal direction respectively. The corresponding results 
of transmission IS21I are depicted in Figure 3(c). 

One can see that basically their behavior are simi- 
lar and the peaks of the resonance transmission moves 
to low frequency as the channel length increases, which 
agrees with the previous Fabry-Perot-like oscillation dis- 
cussed in Figure 2(a). Apparently, the current geometry- 
sensible character of the resonance transmission maxi- 
mums in these CSRR patterned channels didn't give any 
direct verification of the ENZ tunneling theory. In cur- 
rent experiments using the configuration of TEM paral- 
lel plate geometry, the transition section were neglected, 
and the observed transmission peaks were more like to 
be due to Fabry-Perot-like oscillations. More work must 
be done to improve the current experiments to realize 
genuine verification of ENZ tunneling. 

Additional, according to our results(shown in Fig. 
3(c)), the fine structures due to the resonance CSRR in- 
clusion also shifts as vary channel length. The departure 
from the retrieving results from one unit cell is very re- 
markable, suggesting the possible failure of the effective 
medium in CSRR inclusions for finite size and its inher- 
ent dispersion properties. 



III. TUNNELING IN HOLLOW 
RECTANGULAR METALLIC WAVEGUIDE 
OPERATED AT THE CUT-OFF FREQUENCY) 



same as the setup of the tunneling experiment in Rcf[4], 
although with the less size in transverse direction than 
adopted in experiments(about 40 cells) due to finite our 
limited computing resource. The details of the CSRR cell 
structures(shown in the inset of Figure 3(a)) are the same 
as in Ref [j, [j| , these CSRR patterns are perfect metallic 
faces on the surface of FR4 substrate(0.2 mm thick and 
e/e = 4.4) in our setup. The retrieval result of effective 
permittivity for one unit cell was shown in the inset of 
Figure 3(a). The details of the retrieval procedure are 
described in RefjlBIJ. 

In Figure 3(b), totally 5x4 CSRR cells are used (5 
in the propagation direction and 4 in the transverse di- 
rection) to form the effective ENZ medium channel, and 
the simulated transmission IS21I vs frequency is depicted. 



A close analogy has been drawn between the propagat- 
ing mode within a waveguide and the plane wave prop- 
agation in equivalent bulk medium with effective con- 
stitutive parameters 41]. Exploiting the inherent disper- 
sion relation of hollow waveguides, people have realized 
plasma-like effective analogous to these of metamaterials. 

For the fundamental transverse electric (TE10) mode in 
rectangular metallic waveguide filled with non-magnetic 
medium, phase constant is as follows: 

(3 = v'fnw/c) 2 - (ir/w) 2 = lu^/Jm^JJ 

where u> is the angular frequency, /io the magnetic perme- 
ability of vacuum and n the refraction index of medium 
filled in the waveguides, w is the fixed width of the 
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waveguides. So by analogy, the effective permittivity e e f t 
can written as, 

e eff = e (n 2 -c 2 /(4f 2 w 2 )) 

c is the speed of light in air and f is the frequency. The 
cut-off frequency /o = cn/wy/e and the effective per- 
meability remains fj,o, independent of the operating fre- 
quency. 

Thereby the TEio mode can be used to emulate the 
propagation of a TEM wave in an classical lossless plasma 
medium. It has been used in artificial plasma simulation 
at microwave frequencies in Reffllj. In 2002, by inserting 
metamaterials with effective negative permeability inside 
the rectangular waveguide below the cutoff (i.e. e e // < 
0), the equivalent double negative metamaterials have 
been realized. 

In a series of recent work by Edwards et alJfJH], many 
state-of-art demonstrations of nearly perfect tunneling 
have been realized, using the dispersive properties of rect- 
angular metallic waveguides to mimic the response of the 
ENZ medium. For the U-type geometry similar to the 
experimental setup in Ref[6j, our full wave simulation 
using HFSS gives nearly the same results (dashed lines 
shown in Fig. 4). Although the above simple equivalent 
analogy to ENZ medium lacks of convincing theoretical 
evidence. It is interesting that the configuration using 
rectangle waveguides can meet well with the ideal ENZ 
theory so well. The nearly perfect transmission around 
the cut-off frequency can be considered as a nice circum- 
stantial evidence of the equivalence of ENZ medium and 
the waveguide. 

Let's reconsider the origin derivation of their theory 
0, Q on the tunneling through narrow ENZ channel, the 
basic starting point is the z component of H (H z ) holds 
as constant in the ENZ channel. So in another word, 
if the -ff z =const condition can be fulfilled in the narrow 
channel in spite of specific filling, the tunneling effect can 
sill occur without the precondition of ENZ medium. 

In the TEio mode in lossless rectangle metallic waveg- 
uides, the H z of TEio mode is written asfTEj . 

■ a ftw . ,n _ 3 
H z =3A 10 —sm(-x)e Pz 

7T W 

where A\q is constant coefficient and ft is the propaga- 
tion constant. Near the cut-off frequency, ft — > 0, and 
ViJ 2 — > 0, So the H z = cons condition are satisfied au- 
tomatically at its cut-off frequency. So all the effects 
predicted by the tunneling theory 0, Qcan be realized. 
And the current present experiment can be thought as 
the more generic demonstration of the original tunneling 
theory of Silveirinha and Engheta. 

From the discussion in Section II, without the transi- 
tion section (LI and L2), an ENZ channel in a parallel- 
plate geometry can't lead to effective tunneling and the 
Fabry-Perot-like resonance peaks appear at certain fre- 
quencies instead . Is it the case for the hollow waveg- 
uide configuration? To investigate this problem, we per- 
form the corresponding simulation with similar configu- 
ration without lateral transition sections as shown in Fig. 
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FIG. 4: Results of transmission coefficients \S2i \ of tunneling 
channels made of hollow metallic waveguides. The model 
is similar to Figure 1(b) with a = 10mm, w = 20mm, and 
o>ch = 0.5mm. All outer faces are PEC boundaries. Solid 
lines with symbols are simulated I&lI for different channel 
Lengths(Z/o = 15, 20, 30mm) using HFSS. Similar solid lines 
without symbols are corresponding results calculated from a 
simplified equivalent circuit model(shown in Figure 6). The 
dashed line is simulated |Sai| of the channel modelled with 
lateral transition sections(Li = L2 = a c h)- 



1(b), in which two PEC plates parallel to the x-y plane 
at replaced the PMC ones as boundaries at x — and 
w(we take w = 20 mm here). The narrow channel are 
filled with air between two metallic rectangle waveguides 
filled with teflon(e ~ 2eo), supporting the propagation 
of TEio mode above its corresponding cut-off frequency 
f = c/2w^/l « 5.3GHz. 

Our simulation results are depicted in Fig. 4. Near the 
cut-off frequency of the air filled waveguide channel(/o = 
c/2w = 7.5GHz), it is surprising that the |S 2 i| curves 
for different channel length(Lo = 15 — 50 mm) shows 
the peaks of the resonance transmission around the cut- 
off frequency. The maximum of corresponding peaks 
approaches unity, i.e. the nearly perfect tunneling are 
achieved, independent of the geometry of the channel. 
No remarkable difference found with the comparison case 
with air-filled lateral transition section. 

Although there are other transmission peak at higher 
frequencies, their resonance frequency are strongly geo- 
metrically dependent and Further analysis of correspond- 
ing field distribution(not shown in this paper) indicate 
they can attribute to Fabry-Perot resonances. 

To confirm the tunneling nature of transmission peak 
near cut-off frequency, the 3D field distribution of the 
electric field vectors in a typical channcl(Lo = 15mm 
and a c h = 0.5mm) at its resonance frequency are plotted 
in Figure 5. The TEio-mode-distributed field are much 
stronger than those in adjacent teflon filled waveguides, 
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FIG. 5: Calculated field distribution of electric vector 
in the rectangle waveguide tunneling channel with length 
Lo=15mm at the resonance frequency near its cutoff. 




FIG. 6: Simplified circuit model for the metallic rectangle 
waveguide channel. 



suggesting the huge squeezing of electromagnetic energy 
by this air-filled narrow channel. In the propagation di- 
rection, there is little change of the electric field which 
the channel. Thereby the phase difference between two 
end of the channel is near zero, as predicted in the ENZ 
tunneling theory Q. 

Only the narrow air-filled channel itself between two 
waveguides can realize the effective tunneling of the 
incident TEio electromagnetic mode. It seems very 
surprising because of the very different geometric and 
impedance mismatch between air-filled channel and two 
teflon-filled waveguide regions. So it is necessary to con- 
sider the E-plane steps discontinuities in metallic waveg- 
uides, which is an classical problem which can be sim- 
plified and analyzed using equivalent circuits model in 
Ref (T5, fl7j ■ To understand our current results, we mod- 
elled the whole system to be a simplified equivalent cir- 
cuit (shown in Figure 6). 

The fundamental TEio mode is propagating in both 



these waveguides with phase constants as follows, 



Pi = V^Mo^i - (tt/w) 2 

where the wave impedance is Zi = u)jj,o/ (3% respectively 
for each section (i=0 for air-filled channel and 1 for teflon- 
filled waveguides). And an ideal transformer [lj] with the 
ratio I : y/s (s = a c h/a) was used in the model to reflect 
the change of impedance. The numerical calculation have 
been described [17} in detail. The calculated S21 for dif- 
ferent Lq have been shown in Fig. 4. The results agree 
very well with those from the full wave simulation using 
HFSS. The difference becomes slightly apparent only at 
high frequency. 

Around the cut-off frequency, the phase velocity be- 
came infinite as (3 — > and The change of E and H is 
along longitudinal direction. So Lo is negligible. On the 
other hand, when s — ► 0, the change of circuit impedance 
will be dominated by the two transformers (with the ratio 
1 : \fs and 1 : y/s) and relatively the effect of the shunt 
admittance can be neglected. (In our calculation, the 
moderate change on the value of the admittance load, a, 
have very weak effect on the position of resonance trans- 
missions.) The super-coupling channel connected two 
waveguides with the same impedance can be obtained. 
The excellent agreements suggest that such a simple cir- 
cuit model may grasp the nature of tunneling physics in 
metallic rectangular waveguides. Further experimental 
demonstrations of the tunneling in waveguides are ex- 
pected in our future work. 



IV. CONCLUSION 

In summary, we systematically investigated the recent 
experiments about the demonstration of the tunneling 
effect through the narrow ENZ channel. It is shown 
that the ability of the electromagnetic waves penetrating 
into the ENZ medium is very necessary in the parallel- 
plate setup, and the transition section is indispensable 
to demonstrate effective tunneling effect. However for 
the experimental configuration made of hollow metallic 
waveguides, the nearly perfect tunneling can be achieved 
without transition section near its cutoff frequency, which 
is a more generic case and exceed the scope of the original 
ENZ tunneling theory. 
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